Multifractal regime transition in a modified minority game model 
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The search for more realistic modeling of financial time series reveals several stylized facts of real 
markets. In this work we focus on the multifractal properties found in price and index signals. 
Although the usual Minority Game (MG) models do not exhibit multifractality, we study here one 
of its variants that does. We show that the nonsynchronous MG models in the nonergodic phase is 
multifractal and in this sense, together with other stylized facts, constitute a better modeling tool. 
Using the Structure Function (SF) approach we detected the stationary and the scaling range of 
the time series generated by the MG model and, from the linear (nonlinear) behavior of the SF we 
identified the fractal (multifractal) regimes. Finally, using the Wavelet Transform Modulus Maxima 
(WTMM) technique we obtained its multifractal spectrum width for different dynamical regimes. 

PACS numbers: 89.65.Gh, 05.45.Tp, 05.10.-a 



I. INTRODUCTION 

In the last years several agents based models for as- 
set returns have been proposed in the literature p], 
One of the reasons for this interest is that the traditional 
models derived from the geometrical Brownian motion 
do not explain adequately many properties of real mar- 
kets. The Minority Game (MG) and its variants consti- 
tute one of the most promising models [3| due to their 
capability of explaining a wider range of properties found 
in price and index signals. It is usual to characterize eco- 
nomic time series from their empirical properties called 
stylized facts [H, which includes multifractal long range 
correlations. There are several ways to characterize the 
long-range correlations from the real time series and from 
its models. Some of these methods are the autocorre- 
lation functions, power spectral densities (either from 
Fourier or wavelets transforms) and probability distri- 
bution functions. In addition, the fractal and the mul- 
tifractal analysis provide more insights, respectively, on 
the self-similar and self-affine scaling exponents. Here 
we use the Hurst exponent obtained from the Structure 
Function (SF) method 0, and the singularity spectrum 
obtained from the Wavelet Transform Modulus Maxima 
(WTMM) method [1 to determine, respectively, the 
fractal and the multifractal structure of signals generated 
by the Minority Game model. 

The MG models are an oversimplified version of com- 
plex systems. Their peculiarity is that several modifica- 
tions on the original model are possible in order to make 
it more realistic and analytically tractable in the sta- 
tionary state via replica or generating functional meth- 
ods [1, H, [l(| EH- Such methods have led to a deep 
understanding of their macroscopic properties. Unfortu- 
nately, many interesting phenomena occur in the non- 
ergodic phase where the analytical approach fails. Here 
we address a version of MG where the synchronicity of 



time transactions is removed [13 . [H 03 allowing agents 
to trade on different time-scales. Remarkably, for differ- 
ent distribution of frequencies, it has been shown that 
these models essentially preserve many of the statisti- 
cal features of synchronized MG versions such as phase 
transition, volatility clustering, fat tail PDFs and so on. 
However, these features now depend on the distribution 
of the trading frequency. In the present work we will 
focus on the case where the amount of information pro- 
cessed by agents is the same. 

This investigation is part of our continuing quest to 
find models that comply with all known stylized facts 
in financial time evolution [l5|. Thus it is important 
to explore the property of multiscaling in the artificially 
generated data. Many empirical studies have indicated 
the presence of multifractal behavior in real data as a 
stylized fact. Surprisingly, MG models incorporate this 
effect through a more realistic trading rules of agents. 

The paper is organized as follows. In Section 2 we 
describe the model, and then we discuss in Section 3 the 
statistical properties of the return time series. In Section 
4 we explain the multifractal concept and the method 
used to detect it. The results of the multifractal analysis 
are discussed in the Section 5. Finally, in the last section 
we present our conclusions. 



II. MINORITY GAME WITH DIFFERENT 
TIME-SCALES 

The simplest version of the Minority Game [l6[ con- 
sists on a set of N s adaptive agents, also called specula- 
tors. They are endowed with S strategies that map public 
information ji 6 [0, . . . , P — 1] to the decision of buying 
or selling assets, af s (t) — ±1, i — 1 . . . N s , s — 1 . . . S. 
Each strategy is generated at the beginning of the game 
and is kept frozen along the dynamics. Besides, the 
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strategies have a payoff function, Ui >s {t), that describes 
their performance all the time according to the following 
payoff function: 

U itS (t+l) = U i<s (t)-af tS (t)A(t), (1) 

where A(t) — ^ i s , (t) is the excess demand. This rule 
is applied to all strategies independently of their previous 
scores. The public information is a non negative integer 
number drawn at random each round uniformly from in 
the interval [0, P). The decisions are made according to 
af s , (t), where s* is the label for the best ranked strategy 
(s* = max s Ui. s ). In case of a tie, the decision is made 
based on coin tossing. 

An improved model is the grand canonical minority 
game which consists of two distinct groups, namely, pro- 
ducers and speculators. The producers have one strategy, 
their decisions are a function of fi only, and hence they 
are always doing transactions in the market. The spec- 
ulators now have an extra strategy, s — 0, that allows 
them not to trade when their strategies are not prof- 
itable enough. We compute the strategies' performance 
with the following payoff function 

U itS (t + 1) = U it ,(t) - aV a (t)A(t) + e5 , s (2) 

where e is a small real number, which can be positive (in- 
terest rate) or negative (agent's risk aversion measure). 
Thanks to the null strategy, the number of traders fluc- 
tuate along the time evolution. 

Agents in the standard MG models trade on the same 
time-scales and different events also occur with the same 
frequency. We can improve this simplified assumption 
by introducing g groups of N s speculators, each one with 
different time scale tsj, where one group gj play a sin- 
gular minority game or grand canonical minority game 
for j — 1 . . . g. All agents, independently of the group, 
have access to the same information /i and their pay- 
off function are updated virtually all the time. We will 
implement the grand canonical version, with only one 
producer's group with N p members taking action in the 
smallest timescale. Those timescales are chosen to be 
daily, weekly and monthly to resemble the real market. 
The order parameters take into account the total number 
of agents and it is written as a = j^. 

Other authors [H, 0, EH have studied similar mod- 
els. The basic picture of the behavior of the MG, even 
considering different distribution of trading time-scales, 
remains analogous to the synchronized trading frequency 
model. The control parameter a exhibits some depen- 
dence with the distribution frequency parameter. When 
the time-scale difference \tsi~tsj \ increases, for i j, the 
groups become independent and behave as a monochro- 
matic MG. Curiously, the market has the largest fluctua- 
tions for small time-scale difference, i.e., when the groups 
process data that are very close in time. 



III. STATISTICAL ANALYSIS 

The minority game models and their variants were ex- 
tensively explored in several studies. Their capability 
to generate universal statistical properties qualitatively 
similar to the financial market have been the main rea- 
son for this interest. From the statistical point of view, 
the known fact is that the probability distribution of re- 
turns in the minority game time series is Gaussian ex- 
cept around a critical point. In the subcritical region it 
is possible to find certain realizations with fat tails but 
not necessarily following a power law decay. When one al- 
lows interaction between groups composed by the players 
trading in different time scales, the statistical properties 
are qualitatively preserved. 
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FIG. 1: Top: phase diagram for a and H function. Bottom: 
Kurtosis. The critical value of a occurs for the minimum of 
the a curve, around a c — 0.2. The error bars were obtained 
from 40 realizations of different initial conditions. 

Fig. [TJa) presents the variance a 2 and the predictabil- 
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ity H (as a function of a) of the return A(t) defined as: 

fi=0 

where (...) is the temporal average. From this plot we 
can see the phase transition. The critical a c separates the 
ergodic phase from the nonergodic phase. For a > a c 
the time series generated by this model has a normal 
distribution. Fig. QJa) shows the plots of a 2 (upper plot) 
and H (lower plot) against a. Qualitatively these these 
diagrams agree with the synchronous MG model and also 
they help to locate the value of a c . Using the criterion 
that a c is the minimum of a we find that this value is 
near 0.2. In Fig. QJb) we exhibit the kurtosis, denoted 
by k also as a function of a. Firstly in the low a region, 
many realizations have high kurtosis and large error bars 
around a = 0.02, caused by eventual coordination among 
agents. It can be shown that for small values of a the 
time series is leptokurtic. 

The next measure we looked at were the probability 
density functions (PDF) of the fluctuations at different 
scales. These are simply histograms of the height differ- 
ences 

AY r =Y{t + r) -Y(t), (4) 

where Y(t) = J2t'=o A{t') and t is the scale. We will also 
refer to r as the scale of analysis. 

In the Fig. [5] it is shown the PDFs for two values of a, 
for different dynamical regimes. In the bottom panel, for 
a = 0.84, the two groups play independently from each 
other, and the PDFs for all scales are gaussians. In the 
top panel, for a — 0.05, the two groups play interactively 
and the PDFs are fat tailed for small scales. This is a 
sign of intermittency frequently found to be related to 
multifractal processes (l7| . 

IV. STRUCTURE FUNCTION ANALYSIS 

It is usual to characterize the MG dynamics only from 
the statistical analysis point of view as was done in the 
section IIIII In the following we introduce the Structure 
Function (SF) analysis as a preliminary study to the mul- 
tifractal scaling. 

SF have been largely used in the study of turbulence 
[TtJ and are defined in the following way Q : 

S q (T) = (\Y(t + T)-Y(t)\^<x^"\ (5) 

where ( ) denotes the ensemble average. The SF can be 
regarded as a generalization of the correlation functions 
(when q — 2). A signal that is scale invariant and self- 
similar is called fractal when h(q) is the same for all q, 
otherwise it is multifractal [Hj]. Another feature of the 
SF method is its capability to identify nonstationarity 
in the data: for stationary time series the exponent of 
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FIG. 2: Probability distribution functions for two values of 
a. Top: for a — 0.05 the PDFs for small scales are fat 
tailed. Bottom: for a — 0.84 the PDFs for all scales are 
gaussians. Shown as reference, the continuous line, the up- 
side down parabola in both panels are gaussians PDFs with 
zero mean and unit standard deviation. 

the SF is zero, due to the translational invariance of all 
statistics. 

The time series analysis of the standard MG was done 
in [l5j using the unit root test and it was found that the 
MG time series returns are stationary. To extend the 
analysis to the present model, we first verified that the 
SF of the returns, A(i), for q in the range [1, 6.5] were all 
flat for time series 100, 000 points long, indicating that 
the signal is stationary. Although the time series for the 
returns was found stationary in a broad range, ones still 
needs to determine the scaling regime for the cumulative 
sum of the returns, Y(t), since its scaling range does not 
necessarily follows the stationary range. 

This difference between the stationary range and the 
scaling range is seen in the SF of Y(t) for a = 0.05, panel 
a) in the Fig. H These curves are clearly flat for r > 1000. 
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FIG. 3: Top row: structure functions, S q (r), for the signal Y(t), in arbitrary vertical scales, with dashed vertical lines indicating 
the scaling range. From top to bottom, the curves are for q = 1.0, 2.0, ...6.0. Bottom row: the scaling of the structure functions 
as given by Eq.[5] when h(q) is not constant with q, the signal is said multifractal (bottom, left panel), otherwise the signal is 
fractal (bottom, right panel). In these panels, the continous line slope gives h(q = 2), the Hurst exponent. The left column is 
for a — 0.05 and, the right one, for a — 0.84. 



This does not occur for a = 0.84, panel b) in the same 
figure, that seems to have a scaling range broader than 
for a = 0.05. The superior limit for the scaling range is 
in general different for different values of a: we chose the 
value r = 100, indicated by the dashed vertical line in 
the panels a) and b) , as the superior limit that is met for 
all values of a. While the scaling range limit on the right 
is believed to be due information limits, on the left the 
scaling range is limited by a different reason. Close to the 
interaction scale, r = 5 in the present simulations, that 
is below the scaling range, the system does not have yet 
statistical similarity. We will return to this left limit in 
the next section. In short, panels a) and b) of Fig.[3]show 
clearly that the statistical q order moments follow scaling 
laws. Now we need to identify what is the dependence of 
these scaling exponents with q. 

The bottom line of Fig. [3] shows the scaling of the SF 



as given by Eq. \5\ qh(q), against q. For a = 0.05, panel 
c), the dependence is nonlinear with q, indicating that 
the signal is multifractal, but for a = 0.84, panel d), the 
dependence is linear, indicating that the signal is frac- 
tal. For reference, in the bottom panels, the slope of the 
continuous line gives H, the Hurst exponent. Although q 
may be in principle any real number, negative moments 
are difficult to evaluate due to divergence problems and 
will be treated using the WTMM method in the next 
section. 

V. MULTIFRACTAL ANALYSIS 

A given signal can be considered as self-similar with 
scaling exponent h if its statistical properties are invari- 
ant under simultaneous transformation of time t and am- 
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plitude Y(t) [20j, 
t^bt = t'\ 



Y b h Y = Y'- 



(6) 



where & is an arbitrary positive constant and h is a scaling 
exponent given by the equation ([5]). A usual method to 
compute h is based on the SF approach, as shown in the 
section IIVI 

To obtain the full multifractal spectrum we make use 
of the Wavelet Transform Modulus Maxima (WTMM) 
method. The wavelet family used in this paper is the 
n* ft -derivative of the Gaussian density function(DOGn), 
whose wavelet transform has n vanishing moments and 
removes polynomial trends of order n— 1 from the signal. 
Because the scaling properties of the signal are preserved 
by the wavelet transform, it is possible to obtain its mul- 
tifractal spectrum using this method. The number of 
vanishing moments for the wavelet basis, n, is chosen to 
match the order up to (n — 1) of the polynomial trends 
in the signal. 

The wavelet transform of a signal Y(t) is defined as: 

r^(r,M = -i;iW*(— ). ( ? ) 

4=1 

where r > is the scale being analyzed, tp is the mother 
wavelet DOGti and N is the number of points in j direc- 
tion. In this paper we used n = 4 for all analyses. 

The statistical scaling properties of the singular mea- 
sures found in time series can be characterized by the 
singularity spectrum, D(h), of the Holder exponents, h. 
A possible approach to obtain the singularity spectrum 
directly from the time series is the WTMM method 0, 0] • 
The singularity spectrum, D(h), and the Holder expo- 
nents, h, are obtained from the time series by the follow- 
ing equations: 



h(q) = 



lim n T— V %{q;T,bi(T)}\a.\T^[T,bi(T) 
r^o mr ' — ' 



{bi(r)} 



lim Z(q;r) 

t^o inr 



(8) 



D(h) 



where 



lim 

(6i(r)} 

lim Z*(q;T) 

i-^o m r 



T^[q;T,bi(T)] = 



T^[q;r, 6 l (r)]ln \T^[q;r, 6,(r)]| 

(9) 



\T^[t, 6i(r)]| 



E {bl{ r)}\TAr,h(r)}\ q 



(10) 



and the sum is over the set of the WT modulus maxima 
at scale r, {^(r)}. The singularity spectrum, D(h), 
and the Holder exponents, h, are obtained from scaling 
range of the plots of Equations and © , indicated by 
the dashed vertical lines in the Fig. [U Each curve on 



the figure corresponds to different values of q. For large 
and small scales, the scaling regime is broken: to the 
right, it saturates when the system reaches some physical 
limits and, to the left, the system is in the range where 
the agents interacts. The multifractal spectrum is in the 
Fig. El 
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FIG. 4: The functions Z (top panel) and Z* (bottom panel), 
as given by Eq.[S]and Eq.[5] for a = 0.05. The dashed vertical 
lines show the scaling region used to obtain h(q). Each curve 
corresponds to different values of q. 

One way to interpret the multifractal spectrum in a 
physical sense is by comparison with the Hurst exponent 
for known signals, for instance, the fractional Brownian 
motion [20l l2lj . The fractional Brownian motion can be 
classified according to the probabilities of its fluctuations. 
The usual Brownian motion, obtained from a Gaussian 
white noise, has the same probability of having positive 
or negative fluctuations and H — 0.5. A fractional Brow- 
nian motion with H < 0.5 is more likely to have the next 
fluctuation with opposite sign with respect to last one - 
it is said to be antipersistent. Conversely, a fractional 
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FIG. 5: Multifractal spectra for three values of a — 0.01, 0.05 
and 0.84. These values of a are, respectively, for a multifractal 
regime in the non-ergodic phase, a multifractal regime in the 
transition region and a fractal regime in the ergodic phase. 

Brownian motion with H > 0.5 is more likely to have 
the next fluctuation with the same sign as the last one 
- it is said to be persistent. Antipersistent signals have 
more local fluctuations and seem more irregular in small 
scales. Their variance diverges slower with time than the 
variance of persistent signals. Such signals fluctuate on 
larger scales and look smoother. This discussion is done 
in [22[ and a similar, but more detailed interpretation, is 
given in [23|. 

Now we will analyse for the fractal/multifractal prop- 
erties of the models' dynamics as a function of the con- 
trol parameter a. In the ergodic regime, a = 0.84, Fig. 
[5] shows that the multifractal spectrum is narrow, almost 
collapsed to a single point. As already expected, this re- 
gion is clearly fractal, since the time series for the returns 
is Gaussian. Just before the transition, in the subcritical 
regime, the spectrum is wide, a sign of multifractality. 

The multifractal spectrum may be represented by its 
extrema points, i.e., its minima on the left, hi, and on the 
right, h r , as well as the maximum (top), h$. It is worth 
to note that a wide spectrum (difference between h r and 
hi) is a clear evidence for the multifractal character of 
the time series. 

We studied the behavior of the extrema points as a 
function of a to characterize the transition from the frac- 
tal to the multifractal regime, as shown in the Fig. [5] 
This figure shows that the spectrum width increases as 
a decreases, going from a fractal to multifractal regime. 
This transition is smooth and exhibits large fluctuation 
of hi in the nonergodic phase. The point hi, which shows 
the scaling of the fluctuation of the signal, is sensitive to 
the realization of the trajectory. On the other hand, h r 
captures the scaling of the small fluctuation. The point 
ho fluctuates slightly close to 0.5 which is the value ob- 
served for random walks. 



In the nonergodic, sub-critical phase, the system is said 
to represent an efficient market, where no arbitrage is 
possible: this is coherent with the multifractal nature 
of the time series since this dynamics is richer than the 
fractal one, that is typical on the ergodic phase. 
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FIG. 6: Extrema points of the Multifractal spectra, hi, ho and 
h r , for a values from 0.008 to 0.84. The signature of multi- 
fractality is in the fact that the difference h r — hi increases as 
a decreases below a c — 0.2. In the region a > a c , where data 
are Gaussian, this difference should be vanishingly small for 
time series with infinite number of points. The error bars were 
obtained from 20 realizations of different initial conditions. 



VI. CONCLUSION 

We showed that the MG model presents a very rich 
dynamics with anomalous fluctuations that arise due to 
strong correlation similar to the one observed in systems 
driven out of equilibrium. One of the more difficult styl- 
ized fact to obtain from models for financial systems are 
the long range correlations that generate the fractal and 
multifractal properties present in real time series. The 
synchronous MG model does not present multifractality. 
However, only one ingredient added to the model gener- 
ates multifractal region on the space control parameter 
a. This ingredient was the broken of synchronicity. As 
the statistical properties of the MG are preserved, the 
observed multifractal regime belongs to the nonergodic 
phase, where the market is informationally efficient. 

Although the MG time series are not stationary, their 
increments are. Using the SF approach we detected the 
stationary range of the increments and the scaling range 
from the time series. From the linear (non-linear) be- 
havior of the SF we identified the fractal (multifractal) 
regimes as functions of a. To look at the negative values 
of the moments q we used the WTMM and obtained the 
full multifractal spectrum against a. 
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